c
c     size of the domain
c
      integer nx, ny
      parameter ( nx = 20, ny = 20 )
